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Abstract 

ReacProc is a program written in C/CH — h programming language which can be used (1) 
to check out of reactions describing particles interactions against conservation laws and (2) to 
reduce input reaction into some canonical form. Table with particles properties is available 
within ReacProc package. 

1 Introduction 

ReacProc can be used (1) to check out of reactions describing particles interactions against conser- 
vation laws (for example, for electric charge) and (2) to reduce input reaction into some canonical 
form. 

The followings functionalities were realized: 

1. Parsing of reaction. It means reading reaction from a stream and building corresponding dy- 
namical structure of connection lists. Also reaction syntax is validating against some grammar 
(see section 4); 

2. Reducing of grouping brackets (see section 3); 

3. Particle synonyms processing (see section 3); 

4. Sorting of particle and particle groups in reaction (according to an order in a dictionary) (see 
section 3); 

5. Checking conservation laws and printing messages if it's needed (see section 3); 



6. Print reaction to output stream from corresponding data structure (see section 3). 
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We use PPDL notation [Alekhin and et al. 19876, Alekhin and et al. (1987a)] for reaction. 

Reducing to a canonical form means applying 1, 2, 3, 4, 6 functionalities correspondingly. 

ReacProc is written in C/C++ programming language and uses Bison[Donnelly and Stallman 
2002]+Flex[Paxson 1995] parser generator. 

In order to check conservation laws we need some table with particle properties. This file was 
created from publicly available Particle Data Group file and contains following particle properties: 
electric charge, barion number, S-, C-, B-, T-quantum numbers. There was also created file with 
lepton quantum numbers. ReacProc can check conservation laws corresponding to all of these prop- 
erties. There is also file containing a dictionary. 

To run ReacProc one do call corresponding executable with input file. ReacProc reads line by 
line and process strings ending with "\" symbol is considered as line continuation character 
(multi-fines can be written with it). During processing ReacProc creates several files: 

• "rp-accept.txt" contains reactions satisfying tested conservation laws (in "rp-accept-s.txt" re- 
actions are sorted); 

• "rp-reject.txt" contains reactions which do not satisfy tested conservation laws (in "rp-reject- 
s.txt" reactions are sorted); 

• "rp-unknown.txt" contains reactions which cannot be processed by ReacProc; 

• "rp-log.txt" contains debug messages and additional information about processed reactions. 

Another kind of usage ReacProc (and perhaps more powerful) is usage it as a module in another 
system (for example, in specialized databases). So it's important to describe general structure of 
ReacProc. The rest of the paper is devoted to this point. 

As an example of appHcation let's call possible program for building mathematical observables 
from corresponding reactions (which in their order can be used in some calculations). 

ReacProc available here http://sirius.iliep.su/~siver/rp.zip 

2 Data structures 

The most important data structure is called 'f_s'. Here is its definition written in C programming 
language: 

struct f_s -[ 
struct f_s +1; 
struct f_s *r; 
int c ; 
char *n; 

struct gf_s *deca; 
} 

, where the fields means: 'r,'r' - pointers to abstract left and right (or previous and next) 'f_s' ele- 
ments, 'n' - pointer to name of particle, 'c' - counter of particles, 'deca' - pointer to decay description 
of particle. So, this structure is for one type of particle appearing in reaction. But since it can be 
connected with others such structure, any reference to this structure will be reference to the group 
of particles structures. 

'gf_s' data type is defined this way: 

struct gf_s { 
struct gf_s *1; 
struct gf_s *r; 
struct f_s *pg; 
int sign; 
} 
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, where the fields means: 'r,'r' - pointers to abstract left and right (or previous and next) 'gf_s' 
elements, 'pg' - reference to particle group, 'sign' - unary sign which prescribing to 'pg' group. 

Full reaction is represented as following: 

struct reac { 
struct f_s *is; 
struct gf_s *fs; 
} 

, where 'is' field stands for initial state (sequence of particles types, ordinary in amount of two or 
one), 'fs' field stands for final state (sequence of groups of particle sequences). 

3 Description of the basic functions 

Let's describe most important functions and global variables which they use. 

• The parsing function: 

int yyparse(void) ; 

Description: this function is generated by Bison. It calls 'yylex' function which in its order 
reads reaction from 'yyin' stream. After call 'yyparse' reaction will be in 'r' global variable. 

• The canonization function: 
struct reac* canoni(reac* r) ; 

Description: it transforms reaction 'r' into the canonical form. It calls 'syn_proc_gf_s', 
'reduce_gf_s' and 'sort_full_gf_s' functions. 

• The particle synonims processing function: 
f_s* syn_proc_gf _s (f _s*) ; 

Description: it replaces synonim of particle by the key synonim (which actually is PDG name) 
which stands at first place in synonims Hst of the dictionary file; 

• The function for round brackets (named "grouping") reducion: 
gf_s* reduce_gf _s (gf _s*) ; 

Description: it reduces parenthesis. For example, expression '( EH — h MU+ ) X' is equal to 
'E+ X + MU+ X'. Expression with grouping brackets is represented as "generalized" particle 
with 'n' and 'c' fields equaling to "?compo" and correspondingly; 

• The sorting function: 

gf_s* sort_f ull_gf _s (gf _s*) ; 

Description: it sorts according to 'ord_full_gf_s' and at the lowest level - 'ord(char*,char*)' 
function which has two ordering modes: lexicographical (in case a global 'ord_type' = 
'ORD_TRUE_LEX') that is order that every English dictionary has, and so called dictio- 
nary order (in case the global 'ord_type' = 'ORD_DICT') that is order of some particle 
dictionary (ReacProc uses 'dict-syn.txt' file as a particle dictionary by default). These modes 
have advantages in certain situations: while ORD_TRUE_LEX mode is useful when we com- 
pare reactions, with ORD_DICT mode sorting is more logical and physical well founded (see 
example below); 
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• The function checking conservation law: 

int test_reac(reac *r , const char s [] , prop *pp, int kolvo) ; 

Description: it tests reaction 'r' against law naming 's' and property table 'pp' with 'kolvo' 
elements. It return '0' if no error was detected and some signal (such as RP_SIGN_P_UNK) 
otherwise. It uses global variable 'SIGN_conser' which assigns the type of conservation law (0 
for equality, -1 and 1 for '<=' and '>=' equality respectively. 

• The output function: 

void out_to(FILE *f , reac* s) ; 

Description: it prints reaction 's' into stream 'f. It uses 'out_f_s' and 'out_gf_s' functions. 

Let's write a reaction and see how it's represented. 

E+ E > W- < QUARK QUARKBAR > W+ < TAU+ NUTAU + MU+ NUMU + E+ NUE > 

, where brockets denote sub-decays processes. Using the labels (Ax, Bx, Cx, Dx, Ex) instead of 
memory pointers memory structure corresponding to the reaction can be shown by following ('->' 
sign shows value of the variables and '-->' shows value of the references) 
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After canonization in the case of 'ord_type=ORD_TRUE_LEX' the reaction has following form : 
e+ e- --> W+ < e+ nu(e) + mu+ nu(mu) + nu(tau) tau+ > W- < QUARK QUARKBAR > ; 
and in the case of 'ord_type=ORD_DICT': 

e+ e- --> W+ < e+ nu(e) + mu+ nu(mu) + tau+ nu(tau) > W- < QUARK QUARKBAR > ; 

(notice order of 'tau+' and 'nu(tau)'). 

4 Appendix A. The grammar of reactions 

This appendix contains reactions grammar using in ReacProc. Entry point is 'reaction'. It was cut 
out from corresponding Bison *.y file. 
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reaction : iiiitial_state ARROW fiiial_state END 
I error END 

s 

initial_state : beam target 
beam : PARTICLE 

s 

target : /* empty */ 
I PARTICLE target 

9 

particle : PARTICLE 

I PARTICLE LC final state RC 

I LP final state RP 

s 

decay _group : /* Empty */ // d._g. has type f_s 
I particle decay_group 

s 

final_state : decay _group // f .s. has type gf_s 
I decay _group PLUS final state 
I decay _group MINUS final state 
I final state PLUS CC 

Terminal elements (such as ARROW, END and etc) are described below. This Flex code was 
cut out from corresponding *.lex file: 



"-->" return(ARROW) ; 





return (END) ; 


(" 


return (LP) ; 


)" 


return (RP) ; 


<" 


return (LC) ; 


>" 


return (RC) ; 


+" 


return (PLUS) ; 




return (MINUS) ; 


CC 


' return (CC) ; 



[ \t\n\r]+ ; /* Do nothing */ 

[A-Za-z\+\-\:\"\.\*\=\7.\_\(\)0-9\/]+ return (PARTICLE) ; 
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